Origin of Intrinsic Gilbert Damping 
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The damping of magnetization, represented by the rate at which it relaxes to equilibrium, is 
successfully modeled as a phenomenological extension in the Landau-Lifschitz-Gilbert equation. 
This is the damping torque term known as Gilbert damping and its direction is given by the vector 
product of the magnetization and its time derivative. Here we derive the Gilbert term from first 
principles by a non-relativistic expansion of the Dirac equation. We find that this term arises when 
one calculates the time evolution of the spin observable in the presence of the full spin-orbital 
coupling terms, while recognizing the relationship between the curl of the electric field and the time 
varying magnetic induction. 
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The Gilbert damping torque in magnetic systems de- 
scribes the relaxation of magnetization and it was intro- 
duced into the Laudau-Lifschitz equation [TJ [5] for de- 
scribing spin dynamics. Gilbert damping is understood 
to be a non- linear spin relaxation phenomenon and it con- 
trols the rate at which magnetization spins reach equilib- 
rium. The introduction of this term is phenomenological 
in nature [3] and the question of whether it has an in- 
trinsic physical origin has not been fully addressed, in 
the face of rather successful modeling of the relaxation 
dynamics of measured systems. Correlating ferromag- 
netic resonance spectral line-widths |H [5J in magnetic 
thin films with the change in damping has been success- 
ful for confirming the form of the damping term in the 
underlying dynamical equations. The intrinsic origin of 
the damping itself is still an open question. The damping 
constant, a is often reformulated in terms of a relaxation 
time, and the dominant relaxation processes are invoked 
to calculate this, but this approach presupposes preces- 
sional damping torque. 

It has been long thought that intrinsic Gilbert damp- 
ing had its origin in spin-orbital coupling because this 
mechanism does not conserve spin, but it has never been 
derived from a coherent framework. Non-local spin re- 
laxation processes [6] and disorder broadening couple to 
the spin dynamics and can enhance the Gilbert damp- 
ing extrinsically in thin films and heterostructures. This 
type of spin relaxation, which is equivalent to ensemble 
dephasing [7], is modeled as the (S-So)/T2 decay term 
in the dynamical Bloch equation, where Tj is the decay 
time of the ensemble of spins. Crudely speaking, during 
spin relaxation, some spins lag behind the mean mag- 
netization vector and the exchange and magnetostatic 
fields then exert a time dependent torque. Calculations 
on relaxation driven damping of this kind presuppose the 
Gilbert damping term itself which begs the question. 
The inhomogcncous damping term can be written as 
M x d\7 2 ~M./dt which gives rise to non-local effects such 
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as spin wave dissipation [BJ [5]. These non-local theo- 
ries are successful in quantifying the enhancement of the 
Gilbert damping, but do not derive the intrinsic Gilbert 
term itself. There are models [TU] which deal with 
the scattering of electron spins from thermal equilibrium 
in the presence of phonon and spin-orbital interactions 
which is a dynamic interaction and this allows us to de- 
termine the strength of the Gilbert damping for itiner- 
ant ferromagnetic metals, generalizing the Gilbert damp- 
ing response to a tensorial description. Both the s-d 
exchange relaxation models I12j and the Fermi sur- 
face breathing models of Kambersky [HI H3] either pre- 
suppose a Gilbert damping term in the dynamical equa- 
tion or specify a phenomenological Hamiltonian H = - 
l/(7M s )d.dM/dt. While this method is ab initio from 
the point of view of electronic structure, it already as- 
sumes the Gilbert term ansatz. Hankiewicz et al. [14] 
construct the inhomogeneous Gilbert damping by con- 
necting the spin density-spin current conservation law 
with the imaginary part of magnetic susceptibility ten- 
sor and show that both electron-electron and impurity 
scattering can enhance the damping through the trans- 
verse spin conductivity for finite wavelength excitations 
(q 0)- I n previous work [15] . there are derivations 
of the Gilbert constant by comparing the macroscopic 
damping term with the torque-torque correlations in ho- 
mogeneously magnetized electron gases possessing spin 
orbital coupling. For the case of intrinsic, homogeneous 
Gilbert damping, it is thought that in the absence of 
spin-orbital scattering, the damping vanishes. We aim to 
focus on intrinsic, homogeneous damping and its physical 
origin in a first-principles framework and the question as 
to whether spin in a homogeneous time-varying magne- 
tization can undergo Gilbert damping is addressed. 
In this work, we show that Gilbert damping does indeed 
arise from spin-orbital coupling, in the sense that it is 
due to relativistic corrections to the Hamiltonian which 
couple the spin to the electric field and we arrive at the 
Gilbert damping term by first writing down the Dirac 
equation for electrons in magnetic and electric potentials. 
We transform the Hamiltonian in such a way as to write 
it in a basis in which the canonical momentum terms are 
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even powers. This is a standard approach in relativistic 
quantum mechanics and we do this in order to calculate 
the terms which couple the linear momentum to the spin 
in a basis which is diagonal in spin space. This is often 
referred to as a non-relativistic expansion of the Dirac 
equation. This allows us to formulate the contributions 
as a perturbation to an otherwise non-relativistic parti- 
cle. We then wish to calculate the rate equation for the 
spin observable with all of the spin-orbital corrections in 
mind. 

Now, we start with a purely relativistic particle, a Dirac 
particle and we write the Dirac-Pauli Hamiltonian, as 
follows : 



now apply the procedure of eliminating odd powers once 
again : 
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H = e Lb H e~ lb = f3m c z + s + , (5) 

where O is now 0( m l ci ), which can be further elimi- 
nated by applying a third transformation (S = e4 ), 
we arrive at the following Hamiltonian : 



H = ea.fp — e — ) + /3m c 2 + ecf> (1) 
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= O + (3m c 2 + e (2) 

where A and (f> are the magnetic vector potential and the 
electrostatic potential, respectively and 
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We observe immediately that (30 — —0/3. O is the Dirac 

canonical momentum , c and e are the speed of light in 

a vacuum and the electronic charge, respectively. 

We now need to rewrite the Hamiltonian in a basis where 

the odd operators (whose generators are off diagonal in 

the Pauli-Dirac basis : a 1 , 7 l ,7s ..) and even operators 

(whose generators are diagonal in the Pauli-Dirac basis : 

(1, P, X,.. ) are decoupled from one another. 

If we are to find S so that H does not contain odd powers 

of spin operators, we must chose the operator S, in such 

a way as to satisfy the following constraint : 
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In order to satisfy cancelation of the odd terms of O 
to first order, we require S = and this is known 

as the Foldy-Wouthuysen transformation in relativistic 
quantum mechanics and it is treated in some detail in, 
for example, reference [IJj]. We now would like to collect 
all of the terms into the transformed Hamiltonian, and 
this is written as 
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The expression above contains odd powers of the canon- 
ical momentum O, so we redefine the canonical momen- 
tum to encapsulate all of these odd power terms. So we 
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Thus we have the fully Foldy-Wouthuysen transformed 
Hamiltonian : 
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The terms which are present in the above Hamiltonian, 
show us that we have a p 4 kinetic part which is the rela- 
tivistic expansion of the mass of the particle. The terms 
which couple to the spin E are of importance and we see 
that these terms correspond to the Zccman, spin-orbital 
(comprising momentum and electric field curl terms) and 
the Darwin term, respectively. Strictly speaking, the 
presence of the scalar potential 4> breaks the gauge invari- 
ance in the Pauli-Dirac Hamiltonian and a fully gauge in- 
variant theory would require that this contain the gauge- 
free electromagnetic field energy. We omit the term 
4 e 2 " 3 S.(A x E) (which establishes gauge invariance in 
the momentum terms) in this rotated Hamiltonian, as it 
is 0(l/m 2 c 3 ) and we are only interested in calculating 
semiclassical rate equations for fields, which are mani- 
festly gauge-invariant, and not wave functions or energy 
eigenvalues. We can now define the spin dependent cor- 
rections to a non-relativistic Hamiltonian : 
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and Ui are the Pauli matrices. Note that the last 
two terms in Equation [6] encapsulate the entire spin 
orbital coupling in the sense that these terms couple 
the particle's linear momentum to the spin S^. The 
first spin-orbital term in the Hamiltonian is well known 
and give rise to momentum dependent magnetic fields. 
When the ensuing dynamics are calculated for this 
case, it gives rise to spin relaxation terms which are 
linear in spin |17j . Note that, while neither spin-orbital 
term is Hermitian, the two terms taken together are 
Hermitian and so the particles angular momentum 
is a conserved quantity and the total energy lost in 
going from collective spin excitations (spin waves) to 
single particles states via spin-orbital coupling is gained 
by the electromagnetic field. Recognizing the curl of 
the electric field in the last term, we now rewrite this 
the time varying magnetic field as given by Maxwells 
equations as VxE = -^. We now have an explicitly 
time-dependent perturbation on the non-relativistic 
Hamiltonian. We can write the time-varying magnetic 
field seen by the spin (in, for example a magnetic 
material) as = f| • ™ = Mo (l + X™ 1 )^- We now 
have the spin dependent Hamiltonian : 
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We focus our attention on the explicitly time-dependent 
part of the Hamiltonian H s (t) ; 
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In this perturbation scheme, we allow the Hermitian 
components of the Hamiltonian to define the ground sate 
of the system and we treat the explicitly time-dependent 
Hamiltonian (containing the spin orbital terms) as a time 
dependent perturbation. In this way, the rate equation is 
established from a time dependent perturbation expan- 
sion in the quantum Liouvillc description. We now define 
the magnetization observable as M ffA^B TrpS a (t) 

a 

where the summation is taken over the site of the magne- 
tization spin a. We now examine the time dependence of 
this observable by calculating the rate equation according 
to the quantum-Liouville rate equation ; 
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This rate equation governs the time-evolution of the 
magnetization observable as defined above, in the non- 
equilibrium regime. We can write the time derivative of 



the magnetization [IB], as follows ; 
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and we can use the quantum Liouville rate equation as 
defined by Equation [8] to simplify this expression and we 
arrive at the following rate equation : 
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In the case of the time dependent Hamiltonian derived 
in equation [7] we can assume a first order dynamical 
equation of motion given by = 7M x H and calculate 
the time evolution for the magnetization observable : 
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where, in the last two steps, we have used the fol- 
lowing commutation relations for magnetization spins : 
[S^, Sp] = iheijkS^5 a p which implies that the theory pre- 
sented here is that which relates to local dynamics and 
that the origin of the damping is intrinsic. We now rec- 
ognize the last equation as the which describes Gilbert 
damping, as follows : 
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whereby the constant a is defined as follows : 
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The a defined above corresponds with the Gilbert 
damping found in the phenomenological term in the 
Landau-Lifschitz-Gilbcrt equation and x-m is the mag- 
netic susceptibility. In general, the inverse of the suscep- 
tibility can be written in the form [IT?] . 
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where the equilibrium magnetization points along the z- 
axis and oj ex is the excitation frequency associated with 
the internal exchange field. The Sij term in the in- 
verse susceptibility does not contribute to damping mech- 
anisms as it corresponds to the equilibrium response. 
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In the basis (M x ± iM y ,M 2 ), we have the dimensionless 
transverse magnetic susceptibility, as follows : 



Xmi(q,w) = 7Mo 
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The first term in the dimensionless Gilbert coefficient 
(Equation 111 is small (~ 10~ n ) and the higher damp- 



ing rate is controlled by the the inverse of the suscep- 
tibility tensor. For uniformly saturated magnetization, 
the damping is critical and so the system is already at 
equilibrium as far as the Gilbert mechanism is concerned 
(dM/dt = in this scenario). The expression for the 
dimensionless damping constant a in the dc limit (u> =0 
)is: 



In conclusion, we have shown that the Gilbert term, 
heretofore phenomenologically used to describe damping 
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and we have the transverse spin conductivity from the 
following relation (in units whereby h=l) : 
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where rj ls and r^ e are the impurity disorder and electron 
electron-electron scattering times as defined and param- 
eterized in Reference [2] . We calculate the extrinsically 
enhanced Gilbert damping using the following set of pa- 
rameters as defined in the same reference ; number den- 
sity of the electron gas, n=1.4xl0 27 m~ 3 , polarization p, 
equilibrium magnetization M = 7pn/2, equilibrium ex- 
citation frequency ujq = Ep[(l + p) 2 / 3 — (1 — p) 2//3 ] and 
wave- number defined as q = 0.1 , where and 
are the Fermi energy and Fermi wave number, respec- 
tively m* is taken to be the electronic mass. Using these 
quantities, we evaluate a values and these are plotted as 
a function of both polarization and disorder scattering 
rate in Figure [T] 

In general, the inverse susceptibility Xm 1 wm deter- 
mine the strength of the damping in real inhomogeneous 
magnetic systems where spin relaxation takes place, sub- 
bands are populated by spin orbit scattering and spin 
waves and spin currents are emitted. The susceptibil- 
ity term gives the Gilbert damping a tensorial quality, 
agreeing with the analysis in Reference [10] . Further, the 
connection between the magnetization dynamics and the 
electric field curl provides the mechanism for the energy 
loss to the electromagnetic field. The generation of radi- 
ation is caused by the rotational spin motion analog of 
electric charge acceleration and the radiation spin inter- 
action term has the form : 
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FIG. 1: (Color Online) Plot of the dimensionless Gilbert 
damping constant a in the dc limit (u=0), as a function of 
electron spin polarization and disorder scattering rate. 



in magnetization dynamics, is derivable from first prin- 
ciples and its origin lies in spin-orbital coupling. By a 
non-relativistic expansion of the Dirac equation, we show 
that there is a term which contains the curl of the elec- 
tric field. By connecting this term with Maxwells equa- 
tion to give the total time-varying magnetic induction, 
we have found that this damping term can be deduced 
from the rate equation for the spin observable, giving the 
correct vector product form and sign of Gilberts' origi- 
nal phenomenological model. Crucially, the connection 
of the time-varying magnetic induction and the curl of 
the electric field via the Maxwell relation shows that 
the damping of magnetization dynamics is commensu- 
rate with the emission of electromagnetic radiation and 
the radiation-spin interaction is specified from first prin- 
ciples arguments. 



Acknowledgments 



M. C. Hickey is grateful to the Trinity and the uni- 
formity of nature. We thank the U.S. -U.K. Fulbright 
Commission for financial support. The work was sup- 
ported by the ONR (grant no. N00014-09- 1-0177), the 
NSF (grant no. DMR 0504158) and the KIST-MIT pro- 
gram. The authors thank David Cory, Marius Costache 
and Carlos Egues for helpful discussions. 



[1] L. Landau and E. Lifshitz, Phys. Z. Sowiet. Un. 8, 153 
(1935). 

[2] E. M. Lifschitz and L. P. Pitaevskii, Statistical Physics 
Part 2 (Pergamon Press, Oxford, United Kingdom, 
1980). 

[3] T. Gilbert, Magnetics, IEEE Transactions on 40, 3443 

(2004) . 

[4] C. E. Patton, C. H. Wilts, and F. B. Humphrey, J. Appl. 

Phys. 38, 1358 (1967). 
[5] M. Oogane, T. Wakitani, S. Yakata, R. Yilgin, Y. Ando, 

A. Sakuma, and T. Miyazaki, Japanese Journal of Ap- 
plied Physics 45, 3889 (2006). 
[6] Y. Tserkovnyak, A. Brataas, G. E. W. Bauer, and B. I. 

Halperin, Rev. Mod. Phys. 77, 1375 (2005). 
[7] Y. Tserkovnyak, G. A. Fiete, and B. I. Halperin, Appl. 

Phys. Lett. 84, 5234 (2004). 
[8] G. Eilers, M. Liittich, and M. Miinzenberg, Phys. Rev. 

B 74, 054411 (2006). 
[9] V. Kambersky, Canadian Journal of Physics 48, 2906 

(1970). 

[10] D. Stciauf and M. Fahnle, Phys. Rev. B 72 064450 

(2005) . 

[11] S. Zhang and Z. Li, Phys. Rev. Lett. 93, 127204 (2004). 
[12] M. Fahnle, R. Singer, D. Steiauf, and V. P. Antropov, 

Phys. Rev. B 73, 172408 (2006). 
[13] J. Kunes and V. Kambersky, Phys. Rev. B 65, 212411 

(2002). 

[14] E. M. Hankiewicz, G. Vignale, and Y. Tserkovnyak, 

Phys. Rev. B 78, 020404 (2008). 
[15] E. M. Hankiewicz, G. Vignale, and Y. Tserkovnyak, 

Phys. Rev. B 75, 174434 (2007). 
[16] W. Greiner, Relativistic Quantum Mechanics (Springer- 

Verlag, Berlin, Germany, 1987). 
[17] H.-A. Engel, E. I. Rashba, and B. I. Halperin, Phys. Rev. 

Lett. 98, 036602 (2007). 
[18] J. Ho, F. C. Khanna, and B. C. Choi, Phys. Rev. Lett. 

92, 097601 (2004). 
[19] Z. Qian and G. Vignale, Phys. Rev. Lett. 88, 056404 

(2002). 



